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Charged multiplicity distributions and the mean charged multiplicity have been investigated in 
inclusive neutral current deep inelastic ep scattering with the ZEUS detector at HERA using an 
integrated luminosity of 38 .6pb~ 1 . The measurements were performed in the current region of the 
Breit frame, as well as in the current fragmentation region of the hadronic center-of-mass frame. 
The KNO-scaling properties of the data were investigated and the energy dependence of multi- 
plicity distributions was studied using different energy scales. The data are compared to results 
obtained in e + e~ collisions and to previous DIS measurements as well as to leading-logarithm 
parton-shower Monte Carlo predictions. The scaled momentum distributions of charged particles 
in jets have been also measured for dijet photoproduction with the ZEUS detector at HERA using 
an integrated luminosity of 359 pb -1 . The distributions are compared to predictions based on per- 
turbative QCD carried out in the framework of the modified leading-logarithmic approximation 
(MLLA) and assuming local parton-hadron duality (LPHD). The universal MLLA scale, A e s, and 
the LPHD parameter, fc ch , are extracted. 
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1. Introduction 

The HERA ep collider provides a rich field for the study of particle production in a wide range 
of W, the photon-proton centre-of-mass (CMS) energy, and the photon virtuality Q 2 . The data 
presented here were obtained with the ZEUS detector at y/s ~ 300 GeV 2 and concern the study 
of the hadronisation and parton fragmentation processes, phenomena which give deep insight into 
the non-perturbative sector of QCD. Establishing universal features in the properties of the final 
hadronic system in reactions with different initial particles (e + e~, ep, hadron scattering) helps to 
elucidate how the partonic cascades evolve into observed hadrons. In particular, recent results on 
multiplicity distributions in DIS as functions of different energy scales and detailed comparisons 
with e + e~ are discussed below. The formation of hadron jets was also investigated using the 
multiplicity and momentum spectra of charged hadrons in the dijet photoproduction events. The 
measurements verify the validity and consistency of the MLLA approach at energy scale accessed 
at HERA. 

The multipurpose ZEUS detector is described in detail elsewhere [1] 

2. Multiplicity of charged hadrons 

The average multiplicity and multiplicity distributions are being studied intensively in particle 
collisions. In previous studies of DIS events only the virtuality of the exchanged photon, Q, was 
used as the energy scale [2]. A reasonable agreement with e + e~ data was shown except for the re- 
gion of Q below 6 — 8 GeV. Recently, the ZEUS collaboration has performed a detailed study of the 
charged multiplicity in the neutral current deep inelastic scattering (DIS) [3]. These measurements 
of the charged hadron multiplicity are performed in the Breit and in the hadronic centre-of-mass 
(HCM) frames. Due to the restricted detector acceptance only hadrons belonging to the current 
fragmentation regions in both frames were used in the analysis. 
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Figure 1: Comparison of multiplicity distributions Figure 2: Comparison of multiplicity distributions 
in KNO form in bins of 2 • E£ and W (solid markers) in KNO form in bins of M e g- (solid markers) with 
with e + e~ data for one hemisphere (open circles). e + e~ data for the whole event (open circles). 
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ZEUS investigated charged multiplicity 
distributions and mean charged multiplicity 
in terms of different energy scales in order 
to consistently compare ep DIS data with the 
data from e + e~, -vp and |xp scattering. The 
following alternatives to the Q energy scale 
were considered: the energy of the current 
region of the Breit frame E%, the invariant 
mass of the produced particles W, used in 
the current region of HCM, and the invari- 
ant mass of the hadronic system M e g, used 
in both frames. 

The scaling properties of multiplicity 
distributions in a commonly used form, pro- 
posed by Koba-Nielsen-Olsen (KNO) [4], 
were studied in bins of W, 2 • E% and M e fr 
and compared with e + e~ data. The multi- 
plicity distributions in the KNO form from 
ZEUS are shown in Figs. [TJ, ^[ In these plots 
the scaled multiplicity distributions, *P(z) = 
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Figure 3: Mean charged multiplicities as a function 
of M e ff compared to MC predictions (a, b). Compari- 
son of the measurements as functions of M e g with the 
measurements as a function of 2 ■ Eq along with the 
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MC predictions (c). 

{n c h)P(n c h), are plotted as a function of « c h/(«ch)> 
where P(n c h) and (n cn ) are the multiplicity distribution and average multiplicity respectively. Fig. |j] 
shows a comparison of the KNO distributions in bins of 2 • (12 < 2 • E^ < 100 GeV) and in 
bins of W (70 < W < 225 GeV) with measurements in one hemisphere of e + e~, obtained by the 
TASSO collaboration in the energy range 14 < < 44 [5] and by the LEP experiments at 
^/s^ = 91.2GeV [6,7]. There is a remarkable agreement between ep and TASSO data; the LEP 
data differ somewhat from the present measurement in the peak region and at very low n c h/(n c h). 

A comparison of the KNO distributions 
data (both hemispheres) in M e ff 
bins (8 < M e ff < 30 GeV) in the current re- 
gion of HCM and in the current region of the 
Breit frame (8 < M eS < 20 GeV) is shown 
in Fig. ||. There is good agreement between 
the ZEUS and both TASSO [5] and LEP 
data (91.2 < ^/s~ < 209 GeV) [6-8]. The 
mean charged multiplicities were also inves- 
tigated using the different energy scales dis- 
cussed above. In Fig. [3] the mean charged 
multiplicities, measured in the current re- 
gions of the Breit and HCM frames, are pre- 
sented as functions of the invariant mass of 
the corresponding hadronic system and com- 
pared with the MC predictions. All three MC 
models describe the data reasonably well al- 
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Figure 4: Comparison of mean charged multiplicities 
as a function of different energy scales with the data 
from e + e~ and fixed target experiments. 
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though the Herwig prediction is too high in the last bin of M e ff in the current region of the HCM 
(Figs. H] (a), (b)). As is seen in Fig. [| (c), the data in the Breit and HCM frames agree in the region 
of M e ff < lOGeV, while at higher M e ff, (ra c h) ri ses much faster with M e ff in the current region of 
HCM than in the current region of the Breit frame. In Fig. |3](c) the values of 2 • (n C h) as a function 
of 2 -Eft are also plotted. The data follow the same dependence as (n C h) vs. M e ff in the Breit frame 
but differ from those obtained in the current region of the HCM. 

Finally, Fig. ^| shows the comparison of the the mean charged multiplicities in the current 
region of the Breit and HCM frames as a function of 2 • £" B r and W with the data from e + e~ and 
fixed-target experiments. The fixed-target data were scaled by a factor 2 (since they were measured 
in one hemisphere only) and were corrected for the K ( <? and A decays by a factor 1.08, estimated 
using the ARIADNE MC model.The ZEUS measurements show good overall agreement with the 
data from other experiments and exhibit approximately the same dependence on the respective 
energy scale; only the fixed target DIS data deviate at energies above 15GeV. The energy scale 
2 • £g r gives better agreement with e + e~ data at low values of energy than Q. The measurements of 
(n c h) as a function of W agree, within the uncertainties, with the data from e + e~ collisions. 

3. Scaled momentum distributions 

Recently, momentum spectra of charged hadrons in photoproduction were studied in jet frag- 
mentation processes with the ZEUS detector [9]. The results are compared with perturbative 
QCD calculations carried out in the framework of the Modified Leading Log Approximation, 
MLLA, [10] and the hypothesis of Local Parton-Hadron Duality, LPHD [1 l].The MLLA equations 
give an analytical description of the parton shower evolution and an effective scale parameter of 
the QCD calculations, Ags, that is assumed to be universal, i.e. independent of the process consid- 
ered. The LPHD hypothesis predicts that the observed hadron distributions should be related to the 
calculated parton distributions by a normalisation parameter fc ch . Tests of the MLLA predictions 
in conjunction with the LPHD hypothesis permit to expand our understanding of the underlying 
physics of jet fragmentation phenomenon. 

Related studies [12] have been per- 
formed before in e + e~ collisions at LEP 
and PETRA, in DIS ep collisions at HERA, 
(anti)neutrino-nucleon interactions from the 
NOMAD experiment and pp collisions at the 
Tevatron. In the analysis presented here, the 
momentum spectra of charged hadrons are 
studied in dijet photoproduction (yp) events 
from ep collisions. The scaled momentum 
distributions E, = In (1/je p ), where x p is the 
fraction of the jet's momentum carried by the 
charged particle, were measured in restricted 
cones of various opening angles 6 C around 
the jet axis. Jets were reconstructed from energy-flow objects [13] (EFOs) by applying the kj clus- 
ter algorithm [14]. The reconstructed invariant dijet mass was used as an energy scale. It probes 
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Figure 5: yl e ff as a function of fl, where pi denotes the 
characteristic energy scale for each specific process. 
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the range 19 to 38GeV, which spans the energy region between those accessed previously by the 
ZEUS and CDF collaborations. 

To check the validity of the MLLA predictions using the measured t, distributions, two ap- 
proaches were adopted. The first was based on the position of the peak of the £ distributions, ^ pe ak- 
The second was based on a fit of the full shape of the distributions; the limiting spectra, predicted 
by MLLA+LPHD theory [10], were used in the fit in this method. 

In the ^p ea k analysis the values of (§ pea k were extracted from the i§ distributions using a three- 
parameter Gaussian fit. At leading order (LO), the peak position is predicted to be at 
^peak = + VcY — c, where c = 0.29 and Y is a function of the jet energy £j e t and B c (see 
eq.(3) in [9]) and depends also on the parameter A^. Thus the peak position can be directly fit 
to the data, treating yl e ff as a free parameter. The best fit value was found to be yl e ff = 275 ± 
4(stat.)+jj(syst.) MeV for C =O.23. In Fig. |5] the values of A s r are shown as a function of the en- 
ergy scale and compared to the results from different experiments. The data are consistent with the 
prediction that yl e ff is a universal parameter. 



The E, distributions were also fitted 
using the limiting spectrum predicted by 
MLLA. The values of yl e ff extracted from 
these MLLA fits are in reasonable agreement 
with those extracted from the ^ pea k data, al- 
though the values obtained using the MLLA 
fit have larger uncertainties due to sensitivity 
of yl e ff to the choice of the fitting range. The 
values of the LPHD parameters fc ch were ex- 
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Figure 6: K ch as a function of £j e t for the three 6 C 
values. 



tracted also as a function of £j et and B c from 
the fitted limited momentum spectra and are 
shown in Fig. 6. The value of fc ch , measured 
with 8 C = 0.23 and averaged over £j et , was 

CDF collaboration, K^ h = 0.56 ± 0.05(stat.) ±0.09(syst.). The ZEUS data support the predicted 
universality of K^h- 



0.55 ±0.01(stat.)+^(syst.)+^(theo.) and is in good agreement with that reported by 



4. Summary and conclusions 

The charged multiplicity distributions and the mean charged multiplicity have been investi- 
gated in NC DIS ep scattering in terms of different energy scales. Multiplicity distributions in the 
scaling KNO form in the current regions of the Breit and HCM frames exhibit the same behaviour 
as those in one hemisphere of e + e~ collisions when 2 • or W are considered. When energy scale 
M e ff is used, the charged multiplicities exhibit the same KNO-scaling behaviour as those for the 
whole e + e~ event. The energy scales 2 • E£ and W give better agreement with e + e~ data than Q. 

The multiplicity distributions of charged particles in dijet photoproduction events have been 
measured as a function of ^ = ln(l/^ p ). Two methods, the ^ pea k analysis and fit of the ^ distri- 
butions to the MLLA functions, were used to extract the MLLA scale, yl e ff, and LPHD parameter, 
fc ch . The data support the assumption that both parameters are universal. 
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